Formation of the elementary nanocluster in soft matter: an entropic imprint 
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We study the formation of clusters in model protein solutions and colloidal suspensions, using for 
this purpose an accurate treatment in the framework provided by integral equation theories of the 
liquid state. It is well recognized that the existence of a cluster phase in such systems depends on 
the interplay between short-range attraction and long-range electrostatic repulsion in the effective 
interaction among the macromolecular species. Such competing interactions are accounted for in this 
study by a description of molecules in terms of impenetrable spherical particles, interacting through 
a short-range attractive potential followed by a smoother repulsive tail. We show that the entropy 
sensitively signals the appearance of the smallest isotropic cluster (the elementary nanocluster), 
providing an accurate indicator of the threshold temperature whereupon the nature of the system 
changes from a dense gas to a cluster phase. This marked entropic imprint is surprisingly similar 
to that observed in self-assembling systems undergoing gravitational collapse in cosmology. We 
emphasize the broad comprehensiveness of our approach as a predictive tool in soft matter science, 
biophysics and astrophysics. 



Recent small-angle neutron scattering measurements 
and theoretical studies of small volume fraction colloidal 
systems have identified a local peak in the structure fac- 
tor S{q) at a wave vector o much smaller than that of 
the main diffraction peak H, Q . Such a peculiar fea- 
ture differs significantly from what is usually observed in 
"simple" liquids, since the appearance of a local peak 
has been attributed to structural modifications inher- 
ent to the formation of well ordered nanoclusters in the 
fluid phase. In such conditions, molecules in the sys- 
tem are partitioned in a gas of monomers and stable 
nanoclusters of a generic aggregation number. Investi- 
gations in this area have attracted increasing attention 
because they may shed light on aggregation phenomena 
that can precede gelation states, as commonly observed 
in soft matter such as colloidal suspensions and globular 
protein solutions [1, [E, [1, More generally, the con- 
trol of aggregation processes constitutes a crucial step 
in many realms of science and technology, e.g. in col- 
loid and polymer science, in the study of human diseases 
caused by the formation of fibrillar aggregates or protein 
crystals, in the protein structure determination, in nu- 
clear and plasma physics etc. In the outlined scenario, 
probing into the general mechanisms underlying the clus- 
ter formation, and investigating analogies among similar 
behaviors in different fields, are still a challenging task. 
This is particularly true in the field of protein crystal- 
lization 8, 9J, where the definition of systematic proto- 
cols for the optimal control of cluster growth is partly 
hampered by an incomplete understanding of the inter- 
actions among proteins in the mother solution. The de- 
velopment of generic microscopic models, which capture 
the effective interactions in globular protein solutions, 
is made difficult by the complexity and highly special- 
ized properties of such systems (i.e. large anisotropics 
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in the shape and surface charge distributions, sensitive 
dependence on the physico-chemical conditions of the 
environment, etc). In a basic "colloidal" approach, the 
proteins may be described as perfect spheres bearing a 
uniform surface charge, and the effective interaction in 
solution is approximated by a steep (often hard-core) re- 
pulsion followed by a short-range attractive tail of vari- 
ous shapes [13, [Hi 112, E H [H • This model provides 
a good description of the generic phase diagram of both 
colloidal suspensions and proteins [T^ . [T3| . and an ac- 
curate picture of arrested states, including gelation 
and glass transition [l^ , although it cannot fully account 
for the rich behavior of such systems (see e.g. [H, [I3| 
and references therein). Indeed, theoretical studies as 
well as small-angle scattering and confocal microscopy 
experiments in lysozyme solutions and colloid-polymer 
mixtures have revealed that a weaker long-range repul- 
sion must be added to the short-range attraction, in or- 
der to originate aggregates in both systems P, [3, 
The long-range repulsion is attributed to screened elec- 
trostatic interactions, or to the presence of cosolutes in 
solution, whereas short-range forces stem from several 
factors, such as van der Waals interactions, hydrophobic 
effects, depletion mechanisms [sl. H^. [l9l [20| . 



To summarize, experimental evidence and theoretical 
studies on colloids and protein solutions, show which 
minimal requirements a microscopic interaction should 
possess for a nanocluster phase to exist: the latter seems 
to be an appropriate balance between attraction (that 
favors the cluster growth at low temperature) and long- 
range repulsion (that favors the aggregate formation, pre- 
venting a complete phase separation). Therefore, parti- 
cles interacting by means of a two- Yukawa (2Y) potential 
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in the form: 



u(r) 



— cxp 

rja 

H — T- exp 

rja 



-Zi - - 1 



r < a 



r > a 



(1) 



should provide a reasonably flexible and generic model 
to characterize the nanocluster formation in the complex 
fluids at issue from a microscopic viewpoint. More gener- 
ally, 2Y systems exhibit a rich variety of inhomogeneous 
fluid phases and intriguing phase behaviors; as such they 
have been lively studied until recently both theoretically 
and in terms of computer simulation [2l|, [g^, . 2Y 
particles are characterized by the hard-core diameter a, 
the strength and range of the attractive interaction, Ki 
and l/Zi, respectively, and the repulsive strength and 
range, K2 and I/Z2, respectively. For K2 = 0, a. purely 
attractive Yukawa interaction is recovered, while assum- 
ing K2 = Ki = leads to the familiar hard-sphere model. 
The potential in Eq. ^ gives rise to a short-range at- 
traction followed by a long-range repulsion through the 
constraint Ki > K2 and Zi > Z2. 

In this report we provide a theoretical characteriza- 
tion of the transition threshold between the homogeneous 
dense gas and the cluster phase in fluids constituted by 
2Y particles. We focus on the entropy of the fluid across 
the transition temperature, where the microscopic con- 
ditions for the formation of the elementary nanocluster 
(ENC) are first met. Our analysis is carried out in the 
framework provided by the integral equation theories of 
the liquid state (lET). Specifically, in order to determine 
the structural properties of the fluid we have solved the 
Ornstein-Zernike equation coupled with a self-consistent 
closure relation ^24j. The thermodynamic properties may 
then be calculated from the knowledge of both the consis- 
tent correlation functions and the "one particle" bridge 
function [25. . .26J . Entropy is a well recognized indi- 
cator of structural modifications occurring in fluid sys- 
tems [H m, m m, m, iOl. Moreover, given the ab- 
sence of a purely analytical theory, the interpretation of 
experimental results often resorted on the use of integral 
equation theories [l|, [ll|, [H, [13, [Slj. Notwithstanding 
such evidence, as well as the wealth of simulation and 
theoretical studies available on the 2Y fluid, this is the 
first time, as far as we know, that the cluster formation is 
investigated from our viewpoint. As we shall see ahead, 
although usual closure relations of lET yield quite simi- 
lar and accurate predictions for the structural properties, 
this is definitely not the case for thermodynamics. 

The formation of the ENC is analyzed as a function 
of Ki for several values of the ratio K2/K1 of the re- 
pulsion over the attraction strength. This choice of- 
fers the advantage that the interaction potential u{r) 
can serve as a universal function u(r) = e^{r/a), where 
e = {Ki — K2)k^T (with fee the Boltzmann constant and 
T the temperature), being Ki and K2 dimensionless. In 



addition, Eq. ((T]) can be written in a reduced form so that 
(3u{r) = $(r/cr)/r* (where (5 = l/fcer) is a function 
of the reduced temperature T* = [-^^1(1 - 
Hence, increasing Ki while keeping K2/K1 constant, is 
all in all equivalent to lower the effective temperature T* . 
We set Zi — 10, in order to have a 10% ratio of the 
attraction range over the hard-core diameter, as in typi- 
cal protein solutions, Z2 = 0.5, corresponding to a weak 
ionic strength of order 7 = 2.1 mM, and the number den- 
sity p in such a way that the packing fraction (j> = npa^ /6 
is equal to 15%. 

The local structure of the fluid is described in real 
space by the Fourier transform of S{q), i.e. the pair cor- 
relation function g{r), whose peaks indicate the average 
position of successive shells of neighbors around a cen- 
tral particle. The position of second and third shell is 
reported in Fig. lA as a function of Ki, for several fixed 
values of K2/K1. It appears that second neighbors, lo- 
cated at r2{Ki), move progressively towards the first ones 
as the attraction strength Ki increases (or equivalently 
the temperature decreases), until a threshold value 
is reached. Beyond this point, r2{Ki) reaches a plateau 
and remains unchanged within a [iirf';iirf] range of in- 
creasing attraction (indicated in Fig. lA), whereas the 
position of the third shell r^{Ki) further decreases. It 
can be noted from Fig. lA that when Ki is greater than 
Kf, r2{Ki) slightly decreases: the increase of attraction 
and packing effects, due to the closer proximity of third 
neighbors, still induce small local reorganizations around 
the first neighbors. Our theoretical predictions are fully 
supported by standard NVT Monte Carlo simulations: 
the precise agreement between simulation and theory for 
S{q) can be appraised from Fig. IB. The isothermal com- 
pressibility, given by xt = S{q ~ 0)/ pk^T, stays positive 
and finite in all cases under study, indicating that all the 
thermodynamic states lie in a dense vapor region of the 
phase diagram. 

We identify the temperature = \/[K^{l~K2/Ki)], 
corresponding to , as the threshold for the formation 
of the elementary nanocluster. In fact, in the absence of 
any cluster, the "pre-existing form" of the ENC is the 
entity constituted by a central particle surrounded by its 
first and second neighbors; the aggregation process tak- 
ing place at iff may be regarded as the creation of bonds 
between particles of successive shells. Two particles are 
considered as bonded as soon as their average separation 
is less than a little fraction of a and stays almost un- 
changed with increasing Ki , or decreasing T* Q . There- 
fore Ki represents the least attraction (corresponding 
to a low enough temperature), which is sufficient for 
the creation of a bond between first and second neigh- 
bors; then, as shown above, the density of the ENC re- 
mains constant within the range [iff; iff], and increases 
further for higher values of attraction. As visible from 
Fig. lA, iff moves toward lower values as the relative 
height K2/K1 of the repulsive barrier decreases, indicat- 
ing that a progressively higher temperature is sufficient 
to trigger the cluster formation process. It is important 
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FIG. 1: (A) Position of second (circles) and third (crosses) neighbors (in a units) versus K\ for K2/K\ = 0.25, 0.2 and 0.15, 
from bottom to top. For clarity sake, positions for K2/K\ — 0.20 and 0.15 are shifted upward by 0.5 and 1, respectively. Tick 
segments indicate the window [Tff ; iff ], where the position of second neighbors does not change with increasing attraction. 
(B) Comparison between theoretical (solid lines) and simulation (circles) structure factors S{q). Bottom to top: A'l = 1.66, 
2.3 and 2.5 with K2/K1 = 0.25. Results for Ki = 2.3 and 2.5 are shifted upward by 0.5 and 1, respectively. 



to note in Figure IB the progressive enhancement of the 
secondary peak at low-? vectors (that is initially mani- 
fested as an inflection point without a local maximum), 
because it is usually argued that nanoclusters appear as 
soon as the structure factor develops a well defined low-q 
peak. However, this assumption does not appear to be an 
accurate criterion, since the formation of the ENC is well 
underway before the pre-peak itself is clearly displayed. 
As an example, in the case K2/K1 = 0.25, a well de- 
fined low-q peak only begins to appear only at Ki ~ 2.3, 
whereas the ENC is already formed at Ki ~ 1.66. More- 
over, we may recall that some imprecision (error bars) 
in the low-q region may prevent the observation of an 
infiection point in measured structure factors. 

The accuracy of structural results occasions an inves- 
tigation of the formation of the ENC in terms of ther- 
modynamic properties (and specifically entropy) of the 
system. We start from some intuitive considerations. 
The system at issue is composed by a fixed number of 
particles enclosed within a constant volume. In the ab- 
sence of any nanocluster {Ki < Ki), the "pre-existing 
form" of the ENC occupies a given physical volume. At 
the threshold Ki = K^, second neighbors bond with 
the first ones; as a consequence, the volume filled by the 



previous entity is minimum and the space available for 
other particles becomes larger. At this peculiar stage, 
the insertion of any additional particle in the sample be- 
comes a much easier task, in comparison to previous con- 
ditions, and the corresponding energy cost should sud- 
denly decrease. The energy required to insert another 
particle defines exactly the excess chemical potential of 
the fluid, l3ficx- Our theoretical predictions for the lat- 
ter, displayed in Fig. 2A, present the intuitively expected 
behavior: for each K2/K1 value, the excess chemical po- 
tential first increases with attraction, until a jump oc- 
curs precisely at the ENC formation threshold, Xp. Be- 
yond this value, Pficx restarts to increase until a linear 
regime is eventually reached. In Fig. 2A results from 
the Kiselyov and Martynov "one particle" bridge func- 
tion [231 are also shown. The excess chemical potential, 
denoted Pfi^^, increases continuously and monotonically 
and appears completely insensitive to the formation of 
the elementary nanocluster. We point out that, although 
usual closures of the Ornstein-Zernike equation may pro- 
vide refined structural predictions [111 , the observed non- 
responsive behavior of (3lJ^^ exemplifies what one would 
obtain with any usual — but ours — thermodynamic tools 
involved in lET. 
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FIG. 2: (A) Excess chemical potential and (B) entropy per particle (in fee units), as functions of K\ (dotted lines with 
circles) for K^jKx = 0.25, 0.2 and 0.15, from top to bottom. In panel A the vertical lines indicate the iff value. In panel B 
results for K2/K1 = 0.20 and 0.15 are shifted by —2 and —4, respectively; the pair entropy 5(2) (full lines) and the Kiselyov 
and Martynov [23] full entropy Sp^ (crosses) are also shown. 



We may now turn to the notion of order in the fluid 
at issue, and examine the consequences of the above cal- 
culated insertion energies on the total entropy. We know 
that the more the attraction increases (or the tempera- 
ture decreases), the more the system becomes ordered. 
Yet, is it correct to assume that the entropy decreases 
monotonically at the same time? As visible from Fig. 2B, 
the answer is actuafly negative for fluids exhibiting a 
micro-phase separation, such as the 2Y particles at issue: 
as far as Ki is less than K^, bonds are not established, 
and the total entropy S'(t) normally decreases, indicating 
that the system gradually gets ordered; in close proxim- 
ity of K^, the localization of particles around the central 
one progressively increases, enhancing the entropy loss. 
At the threshold K^, the nature of the fluid changes and 
the observed jump in (3^cx has a faithful counterpart in 
the behavior of S'(t) : the loss of available configurational 
space, due to localization, is more than counterbalanced 
from the gain of accessible (physical) volume and there- 
fore the entropy strongly and discontinuously increases. 
We note in Fig. 2B that the total entropy 5*™, calcu- 
lated from /3^™, decreases monotonically all over the 
range of Ki, providing no useful information about the 
ENC formation. Such a monotonic behavior indicates in 
particular that the divergence observed in 5'(t) cannot be 



attributed to the 2Y singularity in Eq. (1) at r = cr. 

We may also analyze the discontinuous behavior of 
5(x) ill the framework of the multiparticle correlations 
expansion of the entropy derived by Nettleton and 
Green[33|. In particular, the pair entropy S(^2)^ given 
by a spatial integral of (functions of) the radial distribu- 
tion function, has been often used as a straightforward 
approximation for the whole entropy (see e.g. [33] and 
references therein). In our case S(^2) decreases monoton- 
ically over Ki (see Fig. 2B): this quantity, "trivially" 
reading the progressive local structuring of the fluid, can- 
not account for the observed jump in 5'(t), that must be 
reflected in the higher-order terms of the expansion. In- 
deed, the behavior of the "residual multiparticle entropy" 
AS*, i.e. the difference between the total entropy and the 
pair entropy (that includes, by definition, the cumulative 
contribution of all correlations involving more than two 
particles), precisely characterizes the structural modifi- 
cations occurring in a fluid of hard spheres as the density 
increases (29J ; more generally, the zero of this function 
has been found to signal the transition threshold between 
less ordered phases and more structured phases in many 
disparate fluids [1^, [H, [H, [s^ . In our case, the plateau 
[iff ; iff], reported in Fig. lA, is characterized by pos- 
itive AS values; 5'(t) crosses S'(2) from the top, leading 
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to AS = 0, systematically at , i.e. at the tempera- 
ture where the ENC density restarts to increase after the 
initial seed has been formed. 

Since the discontinuity discussed above affects the 
thermodynamic potential (the entropy itself in the mi- 
crocanonical ensemble, or the free energy in the canon- 
ical ensemble [stI]), rather than its first derivative, the 
transition may be called of zero order, instead of first 
order. This behavior, at first sight unexpected in col- 
loidal suspensions or protein solutions, is indeed very 
similar to — and supported by — one of the well known 
phase transitions observed in self-assembling systems in 
astrophysics, namely the gravitational collapse. During 
such a zero order phase transition, the onset of aggre- 
gation, as the temperature drops below a critical value, 
corresponds to the starting point of the collapse of the 
second neighbors onto the first shell, after the latter fell 
on the central particle. Statistical physics predicts in 
this context an entropy discontinuity that can be finite 
or infinite (the entropy may diverge to -|-oo) (ssl. I39ll40j . 
The observed similarities are not surprising a posteriori: 
aggregation is a phenomenon common to a wide variety 
of systems, governed by different types of interactions. 
In the present case, whereas self-assembling systems in- 



teract exclusively through attractive forces (possibly in- 
cluding a short-range cutoff), the 2Y particles interact 
through competing forces, where the short-range attrac- 
tion constitutes the dominant contribution. 

We have shown that a neat entropic imprint charac- 
terizes the formation of the elementary nanoclustcr in a 
generic model for effective interactions in complex flu- 
ids, such as globular protein solutions or colloidal sus- 
pensions. Our theoretical predictions involve common 
thermodynamic functions, like entropy or chemical po- 
tential, and are based on an unspccialized description 
of microscopic interactions, but for the presence of com- 
peting short-range attraction and long-range repulsion. 
Nevertheless, we have pointed out that an accurate the- 
oretical approach is essential to reveal a thermodynamic 
signature of the aggregation process, well underway at 
the microscopic level. On the basis of the broad com- 
prehensiveness of our approach, we argue that the ob- 
served entropy behavior should characterize the aggre- 
gation process in disparate systems. The above empha- 
sized close similarities between our study and statistical 
mechanics predictions about zero order phase transitions 
in gravitational collapse, further corroborate this conclu- 
sion. 
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